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ABSTRACT 

The recently discovered apparent dramatic expansion in the effective radii of massive 
elUptical galaxies from z ~ 2toz ~ 0.1 has been interpreted in terms of either galaxy 
mergers or the rapid loss of cold gas due to AGN feedback. In examining the latter case 
we have quantified the extent of the expansion, which is uncertain observationally, in terms 
of the star formation parameters and time of the expulsion of the cold gas. In either case, 
the large global decrease in stellar density should translate into a major drop in the ISM 
density and pressure with cosmic epoch. These cosmological changes are expected to have a 
major influence on the gas accretion mode, which will shift from 'cold' thin disk accretion 
at high redshifts toward 'hot' Bondi fed ADAF accretion at low redshifts. The decline of 
angular momentum inflow would then lead to a spin down of the black hole, for which we 
have calculated more precise time scales; a value of about 0.2 Gyr is typical for a lO^Af0 
central black hole. These results have implications for the different cosmological evolutionary 
patterns found for the luminosity functions of powerful and weak radio galaxies. 

Key words: black hole physics - galaxies: active - galaxies: evolution - galaxies: jets - 
galaxies: ISM - galaxies: interactions 



1 INTRODUCTION 

Powerful radio galaxies (RGs) are known to be hosted by the most 
massive Mq) elliptical galaxies (Matthews, Morgan & 

Schmidt 1964) at all cosmic epochs (e.g., De Brueck et al. 2002; 
Rocca-Volmerange et al. 2004; Seymour et al. 2007; Nesvadba et 
al. 2007). Moreover, the host galaxies are usually found to be ei- 
ther isolated or located in groups as opposed to rich clusters (e.g., 
Longair & Seldner 1979; Best 2004; Hardcastle, Evans & Cros- 
ton 2007). Extensive radio and optical studies have confirmed the 
original inference reached from radio source counts (Longair 1966) 
that, as compared to the present epoch, the space density of power- 
ful RGs and radio-loud quasars was a factor of ~ 10'' higher during 
the 'quasar era' (2 = 2-3) (e.g., Dunlop & Peacock 1990; Willott et 
al. 2001; Grimes et al. 2004). A similarly strong cosmic evolution 
is exhibited by the optically luminous radio-quiet quasars (RQQs) 
(e.g., Hartwick & Schade 1990; Wall et al. 2005) which too reside 
almost exclusively in massive ellipticals (e.g., Falomo et al. 2008). 

In contrast, low power radio galaxies, i.e., those with 
Atsmhz < 1 X 10^^ WHz"^ sr"^ and typically of Fanaroff-Riley 
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(1974) morphology class I,exhibit much weaker cosmic evolution, 
amounting to only a factor of ~ 10 increase in their abundance be- 
tween 2 ~ and 2 ~ 3 (e.g., Jackson & Wall 1999; Willott et 
al. 2001). It has been proposed that the remarkable cosmological 
evolution may be a manifestation of a fundamental change in the 
dominant gas accretion mode powering the AGN since the quasar 
era, which has shifted from a predominantly 'cold' thin disc ac- 
cretion at high redshifts toward 'hot' Bondi type accretion at low 
redshifts (e.g., Hardcastle et al. 2007; cf Cao & Rawlings 2004). 
The hot accretion (dominant at low z) should funnel rather little an- 
gular momentum into the central supermassive black hole (SMBH), 
leading to the possibility of its spinning down. For such RGs, the 
ADAF accretion mode (e.g., Narayan & Yi 1995) is usually invoked 
to explain their observed faint disc emission and low-excitation op- 
tical spectra. Note that while such RGs mostly have an FR I radio 
morphology, this is likely to be dictated mainly by environmental 
factors, especially the ambient density and its gradient (e.g., Gopal- 
Krishna & Wiita 1988, 2000; Snellen & Best 2001; Hardcastle et 
al. 2007; Perucho & Marti 2007; MelUani et al. 2008). 

An observational clue linked to the hypothesis of change in 
the AGN accretion mode with cosmological epoch comes from the 
drastic decrease in the cool gas content of massive ellipticals since 
the quasar era. Observations of CO in nearby massive field ellipti- 
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cals (within 25 Mpc) hiave yielded very low detection rates. Curi- 
ously, the detection rates are higher for their less luminous coun- 
terparts (Sage et al. 2007; Combes et al. 2007; also, Somerville et 
al. 2008). This counter-intuitive result is, however, fully consistent 
with the idea that more massive ellipticals preferentially lost their 
cold gas due to a stronger AGN feedback (e.g., Sazonov et al. 2005; 
Springel et al. 2005; Croton et al. 2006; Bower et al. 2006; Hopkins 
et al. 2008). In contrast, large amounts of molecular gas in high-a: 
massive ellipticals have been found from CO detections (e.g., Pa- 
padopoulos et al. 2000; De Brueck et al. 2005; Klamer et al. 2005) 
as well as PAH detections (Lutz et al. 2008). 

It is interesting that the fraction of sub-millimetre bright RGs 
shows a sharp decline from > 50% at 2: > 2.5 to < 15% atz< 2.5 
(Reuland et al. 2004; also, Greve et al. 2006; Seymour et al. 2007). 
This shows that in RGs the bulk of the expulsion of the cold gas 
and/or its conversion into stars was largely accomplished during 
the quasar era. 

The available data favour the premise that jet-driven gaseous 
outflows were common in RGs during the quasar era and that such 
outflows may have played an important role in the evolution of 
those massive galaxies (e.g., Nesvadba et al. 2008). The mechanical 
feedback of « 5% of the QSO's bolometric power during its peak 
phase can supply the energy needed to induce such intense outflows 
(e.g., Granato et al. 2004; Hopkins et al. 2005). Powerful outflows, 
if produced by nearly-Eddington or even super-Eddington accretion 
onto black holes, have been argued to be fundamentally responsi- 
ble for: the coordinated growth of the masses of the galactic bulge 
and the SMBH (e.g., Silk & Rees 1998; Fabian 1999; King 2003); 
controlling the rate of growth of galaxies (e.g.. Bower et al. 2005; 
Croton et al. 2006; Best et al. 2006); as well as offsetting the cluster 
cooling flows via heating of their ICM (e.g., Binney & Tabor 1995; 
Fabian et al. 2003; Birzan et al. 2004). 

Whereas the expulsion of huge quantities of cold gas from 
massive elliptical galaxies over the past ~10 gigayears may well 
have transformed the nature of their AGN activity, it has recently 
been invoked (Fan et al. 2008) to explain another intriguing result, 
namely the apparent expansion of the effective radii of massive el- 
liptical galaxies by a factor of ~ 3-4 since the quasar era, as found 
in several independent studies (e.g., Ferguson et al. 2004; Trujillo 
et al. 2007; Zirm et al. 2007; van Dokkum et al. 2008; Cimatti et 
al. 2008; Damjanov et al. 2009). These 'superdense' ellipticals are, 
however, found to be exceedingly rare in the local universe (e.g., 
Trujillo et al. 2009). The largest sample of distant galaxies used in 
these studies involved over 800 galaxies with 0.2 < « < 2 and 
masses between ~ 5 x 10^° and ~ 5 x lO^M© (Trujillo et al. 
2007). These ellipticals were strongly argued to be much smaller at 
z ^ 1.5, with the more concentrated (spheroidal) galaxies having 
sizes roughly four times smaller (Trujillo et al. 2007) than those 
found locally in a large sample drawn from the SDSS catalogue 
(Shen et al. 2003). The less concentrated (disk dominated) galaxies 
also showed significant, albeit much more modest, size evolution 
for similarly high masses. 

However a very recent alternative analysis, using a tiny but 
very well observed set of HzRGs, claims that there need be no sig- 
nificant expansion of massive ellipticals and that the multiple ear- 
lier claims of this important effect are in works that did not account 
for the fact that single Sersic profiles are poor matches to the mas- 
sive ellipticals at low z (Hopkins et al. 2009). These authors assert, 
firstly, that they do not find the cores of nearby ellipticals to be 
substantially different in stellar density and size, when compared 
to their high z counterparts. Secondly, the absence of extended en- 
velopes in the existing optical images of high z elUpticals (similar 



to those seen in low 2 ellipticals) could either be a real effect or 
merely an artifact of their low surface brightness (Hopkins et al. 
2009). The former possibility would be consistent with the recent 
work showing that minor mergers would expand the outer portion 
of the elliptical galaxy (e.g., Naab et al. 2009). On the other hand, 
the latter alternative would imply that no substantive change has oc- 
curred in the stellar distribution of the massive ellipticals over the 
past ~10 Gyr that have elapsed since the quasar era. If true, this 
would further accentuate the difficulty in explaining the enormous 
cosmic evolution observed both in the cold gas content of massive 
ellipticals and, even more, in their AGN activity (see above). 

If indeed the evolutionary history of RGs is associated pre- 
dominantly with a transition from a cold disc accretion to a hot 
accretion mode, as mentioned above, a primary requirement would 
be getting rid of most of the cold gas content of the massive ellip- 
tical hosts of the RGs after the quasar era. The changed accretion 
mode would then also result in spinning down of the central en- 
gine. Therefore, a more detailed modelling of both these processes, 
namely the cold gas expulsion and BH spin-down, is of consider- 
able importance and the present paper will address these two issues. 
An improved understanding of the molecular gas expulsion mech- 
anism has now acquired added significance in view of the recent 
controversy over the claimed cosmological expansion of the mas- 
sive ellipticals. 

To further explore these questions, in Sect. 2. 1 we note how 
the possible expansion of the stellar distributions of massive elUp- 
tical galaxies could occur via mergers and in Sect. 2.2 we elaborate 
upon the cold gas expulsion scenario which can lead to the very 
large expansion claimed to have occurred throughout the bulk of 
the galaxies. Different modes of accretion for RGs at different cos- 
mic epochs are considered in Sect. 3. The closely related key issue 
of the spin evolution of the nuclear SMBH in massive ellipticals is 
then discussed in Sect. 4. Given the current debate over the reality 
of the cosmic expansion of massive ellipticals, we shall treat it for 
the present as an unsettled issue while discussing the cosmic evolu- 
tion of the AGN population. However, this uncertainty should not 
affect our main conclusions, which are given in Sect. 5. 



2 GALAXY EXPANSION MECHANISMS 

2.1 Galaxy expansion from galaxy mergers 

Taking the cosmic expansion of massive ellipticals to be a real ef- 
fect, Trujillo et al. (2007) have invoked major 'dry' mergers (those 
in which substantial cold gas is absent). These are expected to dom- 
inate at moderate redshifts over the 'wet' mergers (those with sub- 
stantial amounts of cold gas present) that could be more common at 
z > 2—3 (e.g., van Dokkum 2005). Simulations of major dry merg- 
ers yield few new stars but they do puff up the effective radius of 
the galaxy, roughly as re oc jy^o.^s-i.s^ ^j^j^ stellar mass 

and with the exponent in this relationship declining as the peri- 
centre distance between the conjoining galaxies increases (Boylan- 
Kolchin, Ma & Quataert 2006). Thus, two or three nearly equal 
mass mergers over the past several billion years could account for 
the claimed observed increase in the galaxy size and might also 
explain the prevailing age-uniformity found in local massive ellip- 
ticals (Trujillo et al. 2007). 

But a recent analysis by Bezanson et al. (2009) strongly in- 
dicates that major mergers face greater difficulty in producing the 
reported large expansions and would also probably violate inde- 
pendent constraints on mass growth provided by the evolution of 
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the galaxy mass function. Brighter cluster galaxies (BCGs) have 
been reported to evolve in size with redshift even more quickly than 
most of the early type population; this is most easily understood if 
BCGs grow from many smaller dry mergers (Bemardi 2009). Yet 
another study also indicates that dry mergers have not been impor- 
tant for the evolution of most ellipticals out to 2; ~ 0.7 (Scarlata et 
al. 2007). These different claims may be consistent if results from 
a detailed study of Virgo cluster ellipticals (Kormendy et al. 2009) 
are generally applicable. The brighter ellipticals in Virgo all have 
cuspy cores while the fainter ones do not and these coreless ellip- 
ticals have extra light (above a Sersic profile) in their innermost 
portions. That light could arise from star formation following wet 
mergers, while the more massive cuspy cored galaxies would have 
grown through dry mergers (Kormendy et al. 2009). Both minor 
mergers and dynamical expansion due to mass loss (Sect. 2.2) can 
yield the large increases in galaxy size and also accommodate the 
mass function evolution. They both could also lead to significantly 
more expansion seen in the bulk of the stellar distribution as op- 
posed to just the core (within the central 1 kpc). The "fine tuning" 
in the amount of mass loss needed for that scenario to be effective 
may lead to a preference for the minor mergers as the cause for 
most of the observed expansion (Bezanson et al. 2009). 

A detailed simulation of minor mergers onto a massive ellip- 
tical was recently conducted by Naab et al. (2009). They find that 
the dramatic increase in galaxy size, argued to be present by many 
papers over the past few years, can indeed can be explained in this 
way. Very interestingly, this simulation also shows that the inner 
portion of the stellar distribution remains dominated by the stars in 
the original large elliptical, while the more extended outer reaches 
are mostly filled with the remnants of the smaller galaxies it ab- 
sorbs (Naab et al. 2009). Another recent study of the evolution of 
early-type galaxies indicates that variations in the times at which 
star formation halted for different galaxies can explain about half 
of the observed expansion and dry mergers can explain the remain- 
der of it (van der Wei et al. 2009). Note however, that a recent very 
large sample of 150,000 galaxies matched between SDSS and the 
FIRST radio catalog (Becker et al. 1995), extending to look-back 
times of ~ 2 Gyr, does not show evidence for mergers or other ex- 
ternal environmental factors playing a significant role in triggering 
nuclear activity in either spiral or elliptical galaxies (Reviglio & 
Helfand 2009). 

The merger events of the central black holes are likely to bring 
down the net spin of the new hole due to the random addition of the 
two spins and the spin can be further reduced due to radiation of the 
angular momentimi by gravitational radiation (Hughes & Bland- 
ford 2003). 

2.2 Cold gas expulsion induced galaxy expansion 

An alternative explanation for the apparently observed expansion 
of the effective radii of massive ellipticals by a factor of ~ 3^ 
since the quasar era invokes rapid expulsion of cold gas due to AGN 
feedback in the form of winds and jets (Fan et al. 2008; Sect. 1). 
The presence of large amounts of molecular gas in massive galaxies 
at such high redshifts is indicated by the CO detections (e.g., Pa- 
padopoulos et al. 2000; De Brueck et al. 2005; Klamer et al. 2005; 
Lutz et al. 2008), corroborated by the result that the median redshift 
of sub-millimetre selected galaxies also coincides with the quasar 
era z = 2-3 (e.g.. Chapman et al. 2005). Such massive galaxies are 
typical hosts for bright QSOs (e.g., Falomo et al. 2008), and ob- 
servations of sub-mm QSOs indicate that very substantial amounts 
of cold gas are present in the central regions of their hosts (e.g.. 



Omont et al. 2003, Wang et al. 2007; Lutz et al. 2008). Such sub- 
mm QSOs at redshifts of 2-3 are presumed to be at the transition 
between very rapid star formation and the beginning of powerful 
outflows launched by the Ivuninous central engine which, however, 
may still be shrouded by dust present in the cold gas (e.g., Hopkins 
et al. 2005). The amount of cold gas found in these galaxies at such 
stages is comparable to the stellar mass, of at least the inner several 
kpc (Coppin et al. 2008; Tacconi et al. 2008). 

Direct spectroscopic evidence that is consistent with a very 
powerful outflow in a high z RG (HzRG), with an outflow veloc- 
ity of ~ 1000km s"^ and mass loss rate of ~ 3,OOOM0 yr"^ 
has recently been reported (Prochaska & Hennawi 2009). Likewise, 
good evidence has accumulated for massive outflows of ionized gas 
from HzRGs (z = 2 — 3), possibly due to jet feedback transferring 
around 10% of its power to the outflowing gas (e.g., Nesvadba et al. 
2008; also, Villar-Martin et al. 2007). Also, Tremonti, Moustakas & 
Diamond-Stanec (2007) fotmd strong outflows exceeding 700 km 
s~^ and large mass loss rates in 10 out of 14 galaxies with fading 
starbursts and AGN, albeit with z ^ 0.6. Under suitable conditions, 
even several times weaker outflows would be able to expel the cold 
gas in roughly one Salpeter time (~ 4 x 10^ yr) from the visi- 
ble portion of the galaxy (Fan et al. 2008), though the potential of 
the dark matter halo may well prevent it from being completely lost 
from the galaxy. The feedback of only a small fraction of the QSO's 
total peak power in winds and jets is sufficient to yield such strong 
outflows (e.g., Granato et al. 2004; Hopkins et al. 2005). Recall 
that powerful outflows, if produced by nearly-Eddington or even 
super-Eddington accretion onto black holes have been argued to be 
fundamentally responsible for the coordinated growth of galactic 
bulges and the SMBH mass and for heating of the gas in clusters 
(e.g., SiUc & Rees 1998; King 2003; 2009). 

When the mass loss is rapid compared to the dynamical time 
for the stellar system then it probably causes galaxy expansion 
over the course of several dozen dynamical times, an effect that 
has been long studied for stellar clusters (e.g.. Hills 1980). Fan 
et al. (2008) give an approximate calculation of the galaxy ex- 
pansion for two limiting cases. When the gas ejection timescale, 
Tej <C Tdyn, tiie dynamical timescale, then the initial and fi- 
nal energies are respectively given \sy E = —l3GM^/2R and 
E' = [-/3GM'^/R + I3M'(GM/2R)], where /3 depends on the 
stellar density profile whose shape is assumed to be the same in the 
initial and final configurations, because neither the dispersion nor 
the radius changes during the rapid gas expulsion. 

The ratio of the new to original energies is E'/E = 
{M'/Mf[2 - M/M'] (Biermann & Shapiro 1979). Therefore, 
if M/M' > 2, the system is completely unbound but that is not 
likely to ever occur for massive galaxies. When M/M' < 2 then 
the system will relax to a new equilibrium system where 

a(M, M') = R'/R = (2 - M/M'y\ (1) 

if the new equilibrium system is homologous to the original one. 
Numerical simulations show that when the system is not disrupted 
the new configuration does indeed expand by roughly this factor 
and that it takes about 30-40 initial dynamical times for a new equi- 
librium to be established (e.g., Goodwin & Bastian 2006). In the 
other limit, where Tcj ^ Tdyn, the expansion is through the adia- 
batic orbits of stars, E ^ - ($M + MA$/Tdyn) /2 = 0, which 
implies A$/$ = — Tdyn/Tcj « 0, where $ is the mass averaged 
potential and hence the expansion factor is proportional to the mass 
loss rate (e.g.. Hills 1980) so, R'/R = M/M'. 

Taking a reasonable value of Tdyn w 5 x lO'^yr in the core of a 
massive elliptical the timescale for relaxation should be ~ 2 Gyr so 
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that by 2 ~ 0.8 the compact core galaxies seen at z ~ 2 would have 
settled onto the fundamental plane (Fan et al. 2008). A reasonable 
model would adopt a Sersic profile for the projected stellar density, 
specialized to n = 4, corresponding to a r^^* law (but cf. Hopkins 
et al. 2009), and a Navarro et al. (1997) profile for the total mass. 
With these prescriptions and, for the sake of specificity, adopting 
the evolution model of Granato et al. (2004), the effective radius, 
re, can be found in terms of the virial radius, rvir, the ratio of the 
stellar dispersion to the dark matter dispersion, fa, and the masses 
of the stars cold gas, A/com, and the galaxy halo, Mh, (Fan et 
al. 2008), 



0.34 M* + Mc, 



(2) 



They further argue that reasonable star formation and quasar ac- 
tivity inputs will produce a ratio A/coid/A^* ~ 2/3 for massive 
galaxies nested in halos with Mh > 1O^'^M0. If all of this gas is 
expelled in the fast ejection limit, this would lead to an increase 
in galaxy size by a factor of ~ 3 from Eqn. (1), since M = 
Af* + Afcoid. For less massive galaxies, with Mh < 10^^ Mq, 
the current stellar mass is < 2 x 10^" Af© ; then the nuclear activity 
is probably too weak to eject the gas until a much later epoch, by 
which time Afcoid ^ Mt and then Adcj /A/* where A/cj is the mass 
of the ejected gas, becomes small, too. 

Now we improve upon the estimate of Fan et al. (2008) for 
a{M, M'). To find the factor M/M' we note that initially the mass 
would include the mass in the cold gas Afcoid ; after condensation 
and star formation, part of that mass gets added to A/,. Then the 
ratio of the initial to the final mass, after expulsion of the gas at 
time tx , would be given by 



M ^ M^Mjtx) + M^{t^) 
M' ~ M,{U) 



(3) 



The time begins at the epoch when virialization of the halo is es- 
sentially complete. To make further progress, we invoke the semi- 
analytic model of Mao et al. (2007) using the framework of Granato 
et al. (2004). We have calculated this mass ratio by integrating the 
Eqns A1-A6 of Mao et al. (2007), while making the assumption 
that the expulsion event (which could be either AGN or supernovae 
driven, both of which are accommodated in the formulation) re- 
sults in the entire remaining cold gas exiting from the galaxy. With 
Atburst the length of star formation, the most useful independent 
variable is r = taj/Atburst- Then the mass ratio is found to be 

M _ 1 - exp [-(s7 - l)r5] 

aF ^ ^ 



(s — 1/7) cxp {t5) - 



exp [(57 — l)r(5]/7 ' 



(4) 



where 



&{Mr 



^cond 



^-HM,2)M^^'-\ 



(5) 



with Afi2 the total galactic mass in units of 10^'^ Mq, tcond the gas 
condensation and star formation time, and J-(x) — 1 for x ^ 1 
and J-{x) — for 2; ^ 1. The other quantities in the expres- 
sion for the mass ratio are 7 = 0.7 + O.GA/"]^^'''' [4/(1 + z)] and 
a = tcond/i* ~ 5 for isothermal spheres, with t-^ defined by the 

rate of star formation, Af* ~ — f~~' Af/Af' is primarily a func- 
tion of redshift z, the total mass of the galaxy Af 12, and r, which is 
the time of the gas expulsion expressed in units of the starburst du- 
ration, which is taken to begin at the time the halo became virialised 
(and gas condensation and star formation begin). Plausible values 
have been adopted for the other star formation and condensation 
related parameters (see Mao et al. 2007). 




Figure 1. Values of a, the expansion factor, as functions of the time of the 
cold gas expulsion (measured from the start of star formation in units of 
the starburst duration) for the case of rapid expulsion (r^j <^ Tdy^). The 
four upper curves are for 2 = 5, 3, 2.5, 2, starting from the top. The lowest 
curve is for the case of adiabatic slow mass loss (r^j ^ T-^yn. for the case 
of 2 = 3) where Tgj is the gas exit time scale. For all curves a halo mass of 
M12 = 1 is assumed. 



We would like to point out that Eqn Q, without adding 1, is an 
improvement over Eqn. (5) of Fan et al. (2008) which is an approx- 
imation. The deviation is significant for low r, where additional 
factors in the numerator and denominator of the expression become 
important. As a result, the expansion factor seems to be have been 
underestimated in the model of Mao et al. (2007). In Fig. 1, the ex- 
pansion factors for the case of rapid mass ejection, a(r) are shown 
for a range of the halo virialization redshift, 2 = 5, 3, 2.5, 2; one 
case of slow (adiabatic) mass loss is also shown in Fig.[T] It is seen 
that for a typical value of r = 0.5 — 0.6, and an initial gas-to-stellar 
mass ratio of 2/3, an expansion factor of about 3-4 is expected for 
the rapid expulsion case and less than half as much for the adiabatic 
case. 

We see that while a > 4 is possible it requires rather high 
Afcoid to Af* ratios or very small values of t. Values of r < 0.4 im- 
ply ejection of much of the gas before most of the stars form; how- 
ever, those circumstances are unlikely to yield a massive elliptical 
galaxy. So while this process certainly can puff up the galaxies sub- 
stantially, it requires a large, yet not too large (otherwise the galaxy 
is disrupted), quantity of cold gas to be expelled very rapidly, un- 
derscoring the "fine tuning" problem also noted by Bezanson et al. 
(2009). 

Should the expansion factor be indeed modest with a ^ 1.5 
(Hopkins et al. 2009) there are two alternative possibilities for 
the cold gas depletion. First, as seen from Fig. [T] if the explo- 
sive event is delayed by about 0.1 in r, which is in units of 

8 r 4 1 

Afburst = 6 X 10 J-{Mi2) years, while the starburst 

. 1 2 . 

has been in progress, this would lead to a lower a. There could be a 
variety of reasons for this: either the QSO switches on late or mas- 
sive stars are formed slightly later depending on the details of the 
dumpiness of the gas distribution (we have taken it to be isother- 
mal here). A second possibility envisaged is a steady adiabatic mass 
loss which results in the lower curve of Fig.[T] Either of these two 
situations would be consistent with the negative results of Hopkins 
et al. (2009). 
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3 COSMIC EVOLUTION OF THE ACCRETION MODE 

While the mechanical feedback from the AGN might be responsible 
for the expulsion of the bulk of cold/warm gaseous phase of mas- 
sive ellipticals during the quasar era (Sect. 2), the existing hot gas 
is likely to be largely retained and even replenished by continuing 
stellar mass loss. Thus, independent of the role of cold gas expul- 
sion in causing any expansion of the galaxy, a Bondi type accretion 
of this hot ISM phase onto the central SMBH could become the 
favored mechanism for powering the low-redshift/luminosity RGs 
(mostly of edge-darkened or FR I RGs), particularly because any 
emission lines seen in their optical spectra are usually of the low- 
excitation type (e.g., Hardcastle et al. 2007; Balmaverde, Baldi & 
Capetti 2008; also, Baum et al. 1995). 

In contrast, powerful RGs (FR lis) and RQQs found at high 
redshifts are believed to be powered by a sustained thin disc accre- 
tion of cold gas, which is both capable of spinning up the SMBH 
(e.g., Rees & Volonteri 2007; Berti & Volonteri 2008) and also pro- 
ducing intense disc emission, in addition to ejecting powerful rel- 
ativistic jets whenever conditions near the central engine become 
propitious for jet collimation (e.g., Sikora, Stawarz & Lasota 2007; 
Livio, Pringle & King 2003; Meier 2001). In contrast, the hot ac- 
cretion (dominant at low z) should funnel only little angular mo- 
mentum into the central BH; for such RGs the ADAF accretion 
mode (e.g., Narayan & Yi 1995) can explain their observed faint 
disc emission with low-excitation optical spectra. Evolutionary sce- 
narios along these lines are in vogue (e.g., Sommerville et al. 2008 
and references therein), particularly since such central engines can 
explain the radiatively inefficient "low-excitation RGs" which char- 
acterize the bulk of the FR I population and some less powerful FR 
II sources as well (e.g., Hardcastle et al. 2007; Best et al. 2006; 
also Allen 2006). Only rare instances involving 'wet' mergers at 
low z would push massive ellipticals into the high-excitation AGN 
mode attributable to cold gas accretion. Evidence now exists for the 
hot gas accretion to be the dominant process for FR I RGs across 
three decades in radio power (e.g., Balmaverde, Baldi & Capetti 
2008); that paper reports a nearly linear scaling between the accre- 
tion power and the jet power so that ~ 1% of the rest mass energy 
goes into jet power. 

Thus, it appears that the massive elliptical galaxies that had 
hosted powerful AGN (both radio-loud and radio-quiet) during the 
quasar era will, in most cases, be able to produce only their low- 
luminosity versions at small redshifts. We argue that observational 
support for massive elUpticals undergoing transformation from hot 
to cold accretion mode emerges from the fact that the co-moving 
density of powerful RGs, mostly the edge-brightened FR IPs, peaks 
at (;/) ~ IQ-''-^ Mpc"-' for Fisimhz ^ lO^'^ '^ W Hz^^ sr'^ be- 
tween 2 < z < 3, while the weaker radio sources found at 2 ~ 
(essentially all FR I's) have a comparable density of cj) ~ 10~^'® 
Mpc-3 for PisiMHz ^ 10^* W Hz-i sr'^ (Willott et al. 2001; 
converted to ACDM cosmology by Grimes et al. 2004). Note that 
this radio luminosity is typical of FR I RGs, being nearly an order 
of magnitude below the value for the most powerful FR I sources 
(Best et al. 2005). So, the above scheme for accretion mode change 
would posit that a substantial fraction of the ellipticals hosting the 
FR I RGs at low z are probably descendants of the galaxies that 
had hosted FR n RGs at z > 2, as also hinted by Best et al. (2005). 
This scenario also allows us to vmderstand why the decUne in the 
abundance with cosmic time is so much less steep for low power 
radio RGs than for powerful RGs. 

It is further expected that hot accretion, and thus most FR I 
RGs, would preferentially be associated with more massive galax- 



ies at low 2. This is because the Bondi accretion rate scales ap- 
proximately as MgH, (e.g.. Best et al. 2007). Clearly, hot accretion 
would be favoured by deeper potential wells and this fits with the 
well known result that FR I RGs are usually associated with the 
larger cluster galaxies (e.g., Longair & Seldner 1979; Best 2004; 
Hardcastle et al. 2007) which must also be also richer in hot ISM 
(e.g., Mathews & Brighenti 2003). In the next section we discuss 
some important ramifications of the jet production that arise from 
the switchover to the hot accretion mode at later cosmic epoch. 



4 JET PRODUCTION 

We adopt the flexible paradigm that the jets are either powered to- 
tally by black hole spin, essentially by the Blandford-Znajek (1977; 
BZ) process or also partially from a magnetized accretion disk, as 
in the hybrid model of Meier (1999; 2001). Note that even though 
jets may also form in RQQs, they are liable to be quenched at a 
nascent stage (e.g., Gopal-Krishna, Mangalam & Wiita 2008), and 
hence we would consider only the accretion disk output in the case 
of RQQs. The jet powers in these two cases can be written as (Mc- 
Donald & Thome 1982; Meier 1999) 

{-Cbz = ^ojjrBj^R%j'^c BZ process 
(6) 
£h = ^{B^R^gQ,)^ Hybrid process 

where in the first expression, wf s Qf{^h — Qf)/^h depends 
on the angular velocity of the field, Qf, relative to that of the hole, 
Oif , S± is the component of the magnetic field perpendicular to the 
hole, Rh = r(i)m s (l -h (1 - ff^) GM,/c^ is the radius of 
the event horizon of the BH, where m = GM,/(? , and j = a/m 
is the dimensionless angular momentum of the BH; for the second 
expression, is the azimuthal component of the magnetic field in 
the inner portion of the disk, Rms is the radius of the marginally 
stable orbit, the height of the disk is taken to be J? ~ -Rms, and Q. 
is the angular velocity of the disk at the marginally stable orbit. 

The BZ model is technically derived in the limit of zero ac- 
cretion but the hybrid models, on the other hand, do require at least 
a modest level of accretion, most Ukely in the ADAF mode (e.g., 
Narayan & Yi 1995). As mentioned above, an ADAF type of ac- 
cretion disk close to the SMBH, fueled by essentially Bondi accre- 
tion at larger scales, is now the cormnon premise for the weaker 
"radio mode" for AGN, where jets, but little optical emission, are 
produced (e.g., Somerville et al. 2008). 

Although some general relativistic magnetohydrodynamical 
(GRMHD) simulations indicate that only an essentially BZ type 
of process is likely to yield a relativistic jet while hybrid type mod- 
els do not work well (McKinney 2005), other GRMHD simulations 
do seem to show the flows develop essentially along the line of the 
hybrid models (e.g., Koide et al. 2002; Nishikawa et al. 2005); it is 
fair to say that no consensus has yet emerged on this point. Even 
the BZ process for a massive, rapidly spinning BH does not guar- 
antee that a powerful jet will emerge, since it has been argued that 
successful formation of relativistic jets requires additional collima- 
tion by MHD outflows from accretion disks (Sikora et al. 2007). 
When the accretion rate is close enough to the Eddington rate for 
a "standard" (optically thick but geometrically thin) accretion disc 
to form (e.g., Shakura & Sunyaev 1973), the optical and UV lumi- 
nosity will be high, but powerful jets will only rarely emerge. 

A potentially very interesting outcome of the cold gas expul- 
sion for both the BZ and hybrid models would arise if flux con- 
servation were to be assumed for all the gas that feeds the central 
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engine so that B oc R~ . Then the magnetic field contribution to 
the jet power, which is cx in either case, is enhanced by a factor 
of about a* ~ 100. However, it is more likely that the gas feeding 
the disk, at least initially, is practically within the sole gravitational 
purview of the SMBH and therefore would not partake of the over- 
all gas expulsion. This idea is supported by the observations indi- 
cating much less stellar density evolution in the inner ~ 1 kpc of 
massive ellipticals than in the outer portion of their stellar popu- 
lations (Bezanson et al. 2009). In such an event any effect due to 
magnetic field changes is expected to be much more modest until 
the disk matter is fully accreted by the BH and its spin-down sets 
in. 



4.1 Spin down time scale 

The formulation of the BZ model in Macdonald & Thome (1982) 
leads to the following forms for the jet power, C, and torque, Q 
(details are given in Appendix A), 



£. = Cofr^ij) where £.q = ^^B\g, 



Q = Gojr (j) where Qo 



(7) 
(8) 



where we have dropped the subscript BZ. The geometric factors, 
g and /, are the results of angle averaging over the horizon of the 
magnetic flux and the spin of the field and are model dependent; 
expressions for them are given in Eqns (A6) and (A7). When the 
maximum power is transferred from the BH to the jet, as is typically 
assumed, these factors are of order unity. 

The angular momentum budget and the rotation energy budget 
are given by, respectively 



J 



Joj where Jo = cM,m and 



.1/2 



V2 



where £o ~ M,c 



(9) 

(10) 



For reference we give the numerical values of the various quantities 
used in cgs units: 



Jo = 
Co = 

Go = 



cM.m = 9 X 10®'' M| (g cmVs); 



2 X 10* 



iMi (erg'/s); 



4 X 10*Vs|m| (erg) 



(11) 
(12) 

(13) 



This lets us compute a spin-down time from angular momen- 
tum conservation 



Jo 
Go 



J/ 



7.0 X 10** yrs 



(^(i.,j/)/0.1)] 
BlM9f 



(14) 



where, ji and jf are the initial and final spins, respectively, 
K{ji,jf) is the value of the integral, B4 — B/IO* Gauss and 
Mg = A/./IO^Mq. The spin down time can be found by cal- 
culating the time for the rotational energy to reduce by a factor e; 
see Fig.[2l A detailed evaluation of this factor is given in Appendix 
IBI The time corresponding to e = 1/e with ji = 0.5 is calculated 
to be 0.5 Gyr for A/9 = 1, B4 = 1. 

Another estimate of the spin down time that is more relevant 
observationally is the e-folding time scale of the jet power The jf 
at a time when a reduction by a factor q is reached can be calculated 
from 



dir (ji) =jfr {jf). 




2.0 



Figure 2. The fraction of BH rotational energy, e, lost by the BZ process 
as a function of time is shown for various initial values of the BH spin 
parameter j = 0.2, 0.4, 0.6, 0.8 from right to left. The values assumed for 
other parameters are B/^ = 1, Mg = 1. 




J 

Figure 3. The power of the BZ jet as a function of the BH spin parameter j, 
normalized to the value at j = 1.0. The power goes through a maximum at 
j = before dropping; this arises from the competition between the 

horizon radius and the spin of the BH. 



The above equation leads to a quartic in jf , 

J/ - 2gj/ + = 0, 



(16) 



whose roots are expressible analytically. The positive real root be- 
tween and ji is found to be jf {q, ji), which is fed into Eqn l lB3t 
to calculate the time in terms of tj from Eqn l ll4t . The C{j) goes 
through a maximum at j = \/3/2 (Fig. [3}. The evolution of the 
power as a function of time is given in Fig. ID The time corre- 
sponding to g = l/e,ji = 0.5 is calculated to be 0.5 Gyr for 
A/g — 1,B4 = 1. The spin down time scale is likely to be re- 
duced further when the mass of the hole increases by Bondi accre- 
tion (with negligible net angular momentum) as it spins down. The 

™ ~ (W2)tj ^ ^^^^^ 

(Tj + Ta/2) 

-0.45 Gyr, where Lb is the 



effective timescale is estimated to be r., 



the accretion time scale is Ta 



M.c' 

Le 



(15) 



Eddington luminosity. As a result, Tspin ~ 0.2 Gyr for the typical 
case ji — 0.5, A/9 — 1, /34 — 1. Note that the timescales derived 
above are inversely proportional to the BH mass (Eanll4t. 

The energy loss time-scale is quite similar for the hybrid 
model, although the dependence on BH mass is much more grad- 
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0.0 0.5 1.0 1.5 2.0 2.5 

T (Gyr) 

Figure 4. The power reduction factor, q, of the BZ jet as a function of 
time is shown for various initial values of the BH spin parameter j = 
0.2, 0.4, 0.6, 0.8 from bottom to top. The values assumed for other param- 
eters are B4 = 1, Mg = 1. 

ual. Using the rotational energy from Eqn l llOt above in the hmit 
of low spin, and Eqn (12) for the power in the MHD jet given by 
Meier (1999), we have the time scale 

c w inS -1/10, ,-1/10 . -4/5^-1 

where a is the disk viscosity parameter, m is the accretion rate in 
units of the rate that would produce an Eddington luminosity, and (J 
is a duty cycle parameter expected to be around 10"^ (Meier 1999). 
All quantities are normalized to the power of ten that is indicated by 
the subscript and are chosen to be typical for an ADAF type of flow. 
Thus we find that in either case, for a BH of mass about ~ IO^Mq 
the time-scale for spin-down is a few times lO^yr. Note that in the 
hybrid model the spin-down time scale is extremely insensitive to 
M, (Eqn [17}, while it scales as M."^ for the BZ model (Eqn [14} 
where the angular momentum is directly extracted from the hole 
as opposed to the hybrid case, where it is partly extracted from the 
accretion disk. An important caveat to note is that the spin-down 
time scales derived above (eg. Eqn (14}) depend on the assumption 
that the magnetic field strength is constant during the spin-down. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have focused attention on the transition period 
connecting the quasar era (z — 2 — 3) lo the present era, between 
which the co-moving space density of powerful radio galaxies and 
QSOs has declined by almost three orders of magnitude. As de- 
scribed in the foregoing sections, this "transitional" era is probably 
marked by a changeover in the dominant mode of AGN activity 
from one fuelled by accretion of cold gas onto the SMBH to an- 
other which is dominated by Bondi accretion of hot gas from the 
interstellar medium of the galaxy (Sect. 3, Somerville et al. 2008 
and references therein). Only the cold accretion mode can support 
a strong disc emission and a high spin of the central supermassive 
black hole (SMBH) which can therefore eject powerful jets under 
suitable conditions. An important consequence of this paradigm of 
accretion mode change is that the massive elliptical galaxies that 
were able to host powerful AGNs at high redshifts would usually 
be no longer able to do so at recent cosmic epochs (Sect. 3). 

The postulated basic operational flip of the central engines in 
massive ellipticals since the quasar era underscores the need to 



take a closer look at the physical processes governing the tran- 
sitional era. Clearly, the two dominant processes underlying this 
transformation are the elimination of vast reservoirs of molecular 
gas present in the massive ellipticals at high-z and the consequent 
spinning down of their SMBH. In this paper we have addressed 
both these points quantitatively by presenting an improved theoret- 
ical modelling of the physical mechanisms involved. 

Based on the available evidence we have argued in favour of 
the AGN wind or jet driven expulsion as being the principal mech- 
anism by which massive ellipticals have lost their vast reservoirs of 
molecular gas since the quasar era. Recently, this issue of cold gas 
ejection has attracted a great deal of attention, following the claims 
of an observed ~ 3 — 4 fold increase in the sizes of massive ellip- 
ticals since the quasar era, along with a very recent counter-claim 
(Sect. 1). 

We have obtained a detailed, improved expression for the ex- 
pansion factor, Eqn ([4}, which allows us to distinguish between 
the expulsion histories of galaxies in terms of the star formation 
parameters. The distinction is expressed in our model as the expan- 
sion factor would be lower if more stars form in the core before 
the gas expulsion process begins. Alternatively, a steady adiabatic 
mass loss leads to a modest expansion factor. In the former case, 
the core would be more luminous as more star formation occurs 
and less gas is expelled. More detailed modelling along these lines 
can explore the space of the star formation parameters, the timing 
of gas expulsion, and the cases between adiabatic and rapid mass 
loss, in order to make better connections with observations. 

Irrespective of the potential role in causing the reported dra- 
matic size evolution of massive elliptical galaxies, the cold gas ex- 
pulsion will only have a delayed effect on the output of the central 
engine (accretion disk plus the SMBH). This is because the cold gas 
associated with the nuclear gas cloud that is gravitationally bound 
to the SMBH is expected to be only weakly subject to the AGN 
induced cold gas expulsion. This resulting residual AGN activity 
from the cold accretion era would largely define the transitional era 
and be itself determined by the duration of the spin-down phase 
which the SMBH would trace once the nuclear molecular cloud is 
fully accreted and, consequently, a dramatic fading of the central 
engine has set in. 

We have focused largely on the Blandford-Znajek mechanism 
of extracting rotational energy from the BH in the absence of ac- 
cretion torque. We have derived two time scales using energy loss 
and power reduction as criteria for determining the spin down, 
which turns out to be about 0.5 /{B^Mg) Gyr. If Bondi accretion 
proceeds, it will add to the mass of the BH but little to its angu- 
lar momentum and thus reduce the spin-down time scale to about 
0.2/ {BIMq) Gyr. The evolution of the jet power indicates an in- 
crease before a gradual decline if the initial spin, j > \/3/2, as a 
result of the hole's increasing size. This naturally has implications 
for the evolution of the jet. We plan to expand our work to hybrid 
models in greater detail, and to thus explore disc accretion models 
that explicitly involve angular momentum transport from the hole 
to the disc. 

The transition from cold to hot accretion dominated phase in 
the cosmic evolutionary history of the AGN population is marked 
by a period when the SMBH would continue to be fuelled (and 
spun up) by the accretion of the cold gas located within the sphere 
of influence of the SMBH. This gas is likely to survive the other- 
wise efficient cold gas expulsion due to the intense AGN activity 
during the quasar era. But, even if such nuclear gas cloud has a 
radius as large as 1 pc, it would sustain cold disc accretion phase 
for no more than 10^ years (e.g.. King & Pringle 2007). Moreover, 
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the spin-down of the central engine would also occur in ~ 10® 
years (also, Meier 2001) for reasonable values of mass accretion 
rate and the duty cycle parameter (EgnllVt. Thereafter, i.e., during 
the past 6-7 Gyr, the occurrences of cold accretion dominated AGN 
would become rare, mostly sustained by occasional 'wet' mergers 
of captured galaxies. Thus, the AGN activity over the past several 
gigayears would be increasingly marked by Bondi accretion of hot 
gas powering the central engines and thereby producing mostly FR 
I RGs (Best et al. 2006; Hardcastle et al. 2007). 

To summarize, we have envisaged a scheme characterized by 
the following sequence of events related to the evolution of mas- 
sive elliptical galaxies. At z ~ 3 they contain an abundant supply 
of cold gas that yields both prolific star formation and luminous 
thin disk accretion. The latter usually leads to a fast spinning BH. 
Once a good fraction of the cold gas has been accreted, the result- 
ing fast spinning BHs would become capable of of ejecting power- 
ful jets and forming the luminous radio sources that were so much 
more abundant during the quasar era. But these powerful jets and 
disc winds can easily expel most of the cold gas reservoir, possibly 
resulting in a substantial expansion of the host galaxy. A likely out- 
come of the cold gas expulsion is the eventual dramatic weakening 
of the central engine itself through a spinning down of the SMBH, 
usually on a timescale less than ~ 1 Gyr. 

We thank the anonymous referee for suggestions that im- 
proved the presentation of our results. G-K and RJ.W. thank the 
Indian Institute of Astrophysics for the local hospitality provided 
during their visits and RJ.W. also thanks NCRA for hospitality. 
RJ.W. is supported in part by a subcontract from NSF grant AST 
05-07529 to the University of Washington. 



APPENDIX A: FORMULAE FOR POWER AND TORQUE 



We derive the Eqns l|7Hl3t given in Sect. 14.11 We start from Eqn 
(7.19) in Macdonald & Thorne (1982) for the luminosity through a 
ring at latitutde 9 



Al, = ^ -zu B„Arp, 

4ttc 



(Al) 



where zaj{6) = Rnij) sin 9 is the ring radius, Qnij) is the angular 
velocity of the hole, ^Ip is the angular velocity of the field, Bn = 
B ■ n is the field perpendicular to the hole surface and A>I'(6) is 
the flux through a ring at 9 of arc length RH{j)A9, which is given 
by 

Aip(9) = 2-KTjj{9)RHU)Br,A9. (A2) 
The torque over the same ring is given by 

^^^^ 

Now when the above equations for the power and torque are inte- 
grated over 9 we obtain 



8c ^ 

BlRH{jfnH{j) , 
Tc ^' 

where the geometric factors factors / and g are found to be 



(B-n) 

Bl 

(B-n) 



sm 9 d9 



/2S2f\ 



(A4) 
(A5) 

(A6) 



where B and JIf could be functions of 9 and B^ = 
(B ■ n)^ sin"' 6* d9 is the angle averaged root mean squared 
value of the normal component of the field on the surface of the 
hole. These factors are of order unity; if the maximum power is 
transferred, then Q,f = i^H /2, and both of these integrals are unity. 
Then Eqns ([7][8]l follow after substituting for RhU) and ^IhU), 
which are respectively the black hole's radius and angular velocity 
and are given by 

RhU) = r(j)m = + /)l/^)GM./c^ 



(A8) 
(A9) 



Eqns (9) and (10) then follow as the fundamental equations for the 
black hole angular momentum and rotational energy in terms of 
the irreducible mass of the hole. The coefficients for these physical 
quantities and the luminosity are easily found to be Eqns dl mi3t . 



APPENDIX B: CALCULATION OF THE SPIN DOWN 
FACTOR 

We evaluate the energy radiated when the spin reduces from ji 
jf to be 



AS = CAt = Co 



(Bl) 



As a fraction of the rotational energy budget, this is found to be 



AS 
S 



J 



J/ 



dj 



The integral k defined in Eqn ( 114b is calculated to be 



24m3 



(B2) 



(B3) 



where Wi = r{ji) and Wf = r{jf). Using similar techniques the 
integral 



3f 



(B4) 



Now from Eqn l lB2b . the value of Wf for e = 1/e and ji — 0.5 
is solved for and fed into Eqn ( |B3t . Then using Eqn il4\ the time 
to spin down is found to be ~ 0.5 Gyr. 
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